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^ . ABSTRACT 

' We calculate the gluon dissociation cross-section in an anisotropic quark gluon plasma ex- 

pected to be formed in relativistic nucleus-nucleus collisions. It is shown that the thermally 
weighted cross-section of gluon dissociation undergoes modification in anisotropic plasma 
affecting the J/ip survival probability. The dependence of the cross section on the direc- 
tion of propagation of the charmonium with respect to the anisotropy axis is presented. 
Survival probability of J/ip in two different space time models of anisotropic quark gluon 
plasma (AQGP) has been calculated. It is shown that depending upon the initial conditions 
r^i ! (corresponding to RHIC energies), the survival probability in AQGP differs from that in 
isotropic QGP both in the central as well as forward rapidity regions. For initial conditions 
^ I relevant for LHC energies, marginal difference between the the two space time models has 
OO ' been observed with a given initial conditions. 

vn 
in 

: 1 Introduction 

Ever since the possibility of creating quark gluon plasma (QGP) in relativistic heavy ion 
collision was envisaged, numerous signals were proposed to probe the properties of such an 
^ I exotic state of matter. In this context Satz and Matsui [1] had suggested that the production 
of heavy quark resonances {J/ip) will be suppressed as a result of colour Debye screening in a 
hot and dense system of quarks, anti-quarks and gluons. This suppression could be detected 
experimentally through the dileptonic decay mode of these resonances. ALICE dimuon 
spectrometer [2] is dedicated to look for this type of signal. However, it is a daunting task 
to disentangle the contributions of the heavy quarkonium states to muon spectrum due to 
the background from several other sources, e.g. Drell-Yan, semileptonic decay of open heavy 
flavoured mesons {DD, BB) etc. Low energy muons from kaons and pions also constitute a 
large background. 

In a QGP the much harder gluons can easily break up a J /if) contrary to the case of 
hadronic system. In equilibrating plasma the gluons have much harder momentum sufficient 
to dissociate the charmonium. Such a study has been performed in Ref. [3] quite some time 
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ago. Ever since the physical picture of quarkonium dissociation has undergone shght refine- 
ment during the couple of years. First of all, most of the existing calculations on various 
observables assume, from the very beginning, that the plasma is isotropic which may not 
necessarily be true as we shall argue in the following. Moreover, many properties of the QGP 
are still poorly understood. The measurement of elliptic flow parameter and its theoretical 
explanation suggest that the matter quickly comes into thermal equilibrium (with Ttherm < 1 
fm/c, where Ttherm is the time of thcrmalization) [4]. As for example, one of the major diffi- 
culty is to measure the thermalization {jtheTm) and isotropization (xjso) time of the QGP. On 
the one hand, the success of ideal hydrodynamic fits to experimental data [4] implies rapid 
thermalization of the bulk matter created at RHIC. On the contrary, perturbative estima- 
tion suggests relatively slower thermalization of QGP [5]. However, recent hydrodynamical 
studies [6] have shown that due to the poor knowledge of the initial conditions there is a 
sizable amount of uncertainty in the estimate of thermalization or isotropization time. It is 
suggested that (momentum) anisotropy driven plasma instabilities may speed up the process 
of isotropization [7], in that case one is allowed to use hydrodynamics for the evolution of 
the matter. However, instability-driven isotropization is not yet proved at RHIC and LHC 
energies. 

In absence of a theoretical proof favoring the rapid thermalization and the uncertainties 
in the hydrodynamical fits of experimental data, it is very hard to assume hydrodynamical 
behavior of the system from the very beginning. The rapid expansion of the matter along 
the beam direction causes faster cooling in the longitudinal direction than in the transverse 
direction [5]. As a result, the system becomes anisotropic with {pt^) « {pt^) in the local 
rest frame. At some later time when the effect of parton interaction rate overcomes the 
plasma expansion rate, the system returns to the isotropic state again and remains isotropic 
for the rest of the period. Therefore, it has been suggested to look for some observables which 
are sensitive to the early time after the collision. The effects of pre-equilibrium momentum 
anisotropy on various observables have been studied quite extensively over the past few years. 
Heavy quark energy loss and momentum broadening in anisotropic QGP have been studied 
in Refs. [8, 9]. Effects of anisotropy on photon and dilepton yields have been investigated 
rigorously in Ref. [10, 11, 12, 13, 14]. The effect of initial state momentum anisotropy on the 
radiative energy loss has been demonstrated in Ref. [15]. Recently, the authors in Ref. [16] 
calculated the nuclear modification factor for light hadrons assuming an anisotropic QGP and 
showed how the isotropization time can be extracted by comparing with the experimental 
data. Most importantly, the heavy quark potential has been calculated in Ref. [17] and 
the solutions of Schrodinger equations have been obtained in Ref. [18] corresponding to 
anisotropic system. 

It is to be noted that the calculations of J/ip dissociation cross-section in Ref. [3] have 
been performed in an equilibrating plasma and it is found that the survival probability 
increases in such system. We, in the present work, shall extend the above work assuming 
initial state momentum space anisotropy. 

The plan of the paper is the following. In section 2 we briefly recall the necessary 
ingredients to calculate the thermally weighted gluon dissociation cross section in anisotropic 
media. Then we discuss how this can be implemented to calculate the survival probability 
of J / ip along with space-time models for the anisotropic media. Section 3 will be devoted to 
discuss the results. Finally, we conclude in section 4. 
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2 Formalism 



2.1 The thermal-averaged Gluon-J/^ dissociation cross section 

Peskin and Bhanot first calculated the quarkonium-hadron interaction cross section using 
operator product expansion [19]. Similar result was obtained using the QCD factorization 
theorem in Ref. [20]. Same formalism allows to express the hadron-J/ip inelastic cross section 
in terms of the convolution of the inelastic gluon- J/t/^ dissociation cross section with the gluon 
distribution inside the hadron.The perturbative prediction for ghion-J/ip dissociation cross 
section is [21] 



where g° the energy of the gluon in the stationary J/ip frame; eo is the binding energy 
of the J/ip where go > ^o- Os is the coupling constant and itlq is charm quark mass. A 
few comments about the binding energy of quarkonium states is in order here. It is to be 
noted that We have used the constant binding energy of the J/t/j in AQGP at finite tem- 
perature. However, using the real and imaginary part of heavy quark potential (calculated 
in anisotropic QGP) in Schrodingcr equation the authors of Ref. [22] have shown that the 
binding energy of quarkonium states strongly depends on the anisotropy parameter as well 
as on the hard momentum scale. This observation might have important consequences on 
the gluon dissociation cross-section and hence on the survival probability, 
we assume that the J/ip moves with four-momentum P given by 

P = (Mr cosh y, 0, Pt, Mt sinh y) (2) 

where Mt — ^Mj^^ + P^ is the J/ip transverse mass and y is the rapidity of the J/ip. A 
gluon with a four-momentum K — (A;°, k) in the rest frame of the parton gas has energy 

g° = K.u in the rest frame of the J/ip. Now to calculate the velocity averaged cross section 
in anisotropic media we note that the anisotropicity enters through the distribution function 
[12, 14], 

where Phard is the hard momentum scale, n is the direction of anisotropy which is along the 
beam axis and the parameter ^ is the anisotropy parameter(— 1 < ^ < oo). Phaid is related 
to the average momentum in the partonic distribution function. In isotropic case, ^ = and 
Phard Can be identified with the temperature. In such case the gluon- J /ip dissociation cross 
section becomes [3], 

where I'^i is the velocity between J/ip and a gluon where. 
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Change of variables {K •<-> Q) can be obtained by using Lorentz transformations: 

— [q^E + g(sin 9p sin 9g sin (pg + cos 9p cos 9q)) /Mj/^ (6) 
qE 

^ ^ q+ ipi — [(gMrCOshy - Mj/^)(sin6lpSin6lgSin05 + cos6lpCOs6lg) + |P|]vj/^(7) 

where vj/^ = P/E, P = {E, 0, |P| sin6'p, |P| cos^^p) and q = (g sin 9g cos0g, gsin 9q sin (f)q, gcos ^g) 
In the rest frame of J/iJj, numerator of the Eq. (4) can be written as 

jd'q^a{q')f{k',i,pnara), (8) 
while, the denominator of Eq. (4) can be written as [8] 

j d'kf{k\i,p^^^) = j ci3A;/i,o(^k2 + e(k-n)2,phard) = ^^=L=87rC(3)pLd (9) 

where C(3) is the Ricmann zeta function. The maximum value of the gluon J /ip dissociation 
cross section [21] is about 3 mh in the range 0.7< q^ <1.7 GeV. Therefore high-momentum 
gluons do not see the large object and simply passes through it. On the other hand, the 
low-momentum gluons cannot resolve the compact object and cannot raise the constituents 
to the continuum. 



2.2 Survival probability of J/^ in an anisotropic media 

To calculate the survival probability of J in an anisotropic plasma we consider only lon- 
gitudinal expansion of the matter. With the velocity averaged dissociation cross sections, 
the survival probability of the J/ ip in the deconfined quark-gluon plasma is of the following 
form [3], 

^^^^^ Id^r{R\-r^) 

The upper integration limit Tmax = rnin{T^,Tc) and Tj is the QGP formation time, ^^(t) is 
the gluon density at a given time r. Now the J/%1) will travel a distance in the transverse 
direction with velocity Vj/y,: 



d = —r cos 4> + \J R\ — r'^{l — cos'^cl)) (11) 

Here cos0 = vjj^ ■ r. The time interval = MTd/Pr is the time before J/ip escapes from 
a gluon gas of transverse extension R^. In the case of anisotropic QGP Tc is determined by 
the condition : Phardi^ = Tc) = [23], where Tc ~ 170 — 200 MeV. 

The modifications in Eq.(lO) in anisotropic media come from the gluon density, ng(r) and 



the velocity weighted cross-section. The former is given by n3(r) = 16({3)pl^^^^{T)/[7l'^^Jl + ^{t)] 
and the latter has been discussed in the previous section. The time evolution of ^ and Phard 
is described in the next section. 
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Figure 1: (Color online) The thermal- averaged ghion-J/ip dissociation cross section as func- 
tion of the hard momentum scale at central rapidity {6p = vr/2) for ^ = {0,1,3,5}. (a) 
corresponds to Pt = and (b) is for Pt = 8 GeV. 



2.3 Space time evolution 
2.3.1 Model I 

For an expanding plasma the anisotropy parameter ^ and the hard momentum scale Phard 
(appearing in Eq.(lO)) are time dependent. Thus to calculate S{Pt) one needs to know the 
time dependence of Phard and ^. To obtain the time evolution of the parameters we shall 
follow the work of Ref. [23, 24] and evaluate S{Pt) of the J/ip from the first few fermi of the 
plasma evolution. 

The time dependence of relevant quantities is given by [23], 



Phard(l") 



r 

n 



5(1-A{r)) 



1, 



(12) 



where. 



W(r) 
U 

nix) 



Uir) 



ison5 



U{n) 
1 



3A(r)/4 



l-5(l^A(r))/2 



-[l/{x + 1) + tan"^ ^/Vx] 



(13) 



where the exponent 5 = 2 corresponds to free-streaming pre-equilibrium momentum space 
anisotropy and 5 = corresponds to thermal equilibrium, Tj is the initial temperature of the 
plasma and is the velocity of sound. A(r) = i(tanh[7(r — riso)/Tiso] + 1) is the smeared 
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Figure 2: (Color online) Direction dependence of thermal averaged gluon-J/?/' dissociation 
cross section for two values of Pt and at two different rapidities. 



step function introduced to take into account the smooth transition from non-zero value of 
5 to (5 = at r = riso [14] with 7 being the transition width. 

For isotropic case, we have phard = T, Tiso = Tj so that A = 1, U{t) = Ti/r, and W(rj) = 1. 
By using = 1/3 we recover the Bjorken cooling law [25]. As the colliding nuclei do have 
a transverse density profile, we assume that the initial temperature profile is given by [26] 

T,(r) = T, [2(l-rV/?^)]'^' (14) 
Using Eqs.(12) and (14) we obtain the profile of the hard momentum scale as 

Phard(r,r) = T, [2 (1 -rVi?i)]'^'w'='(r) (15) 

We use two sets of initial conditions for RHIC energies. The first set, henceforth referred 
to as Set I, corresponds to that used in Ref. [3] , i.e., Tj = 550 MeV and Xj = 0.7 fm/c. 
For other set (set II) the initial temperature (time) has been calculated using the measured 
multiplicities at RHIC energies [16] and is given by Tj = 440 MeV corresponding to Tj = 0.15 
fm/c. At LHC energies we use the initial conditions: Tj = 820 MeV and = 0.5 fm/c. 

2.3.2 Model II 

The other alternative scenario of time dependence for C, and phard in highly anisotropic system 
has been described in [27] taking the first two moments of Boltzmann equation which reads 
in (0+l)-dimension as 

df{t,z,p) df{t,z,p) 
E +Pz = -C[/(t,z,p)] (16) 
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Figure 3: (Color online) The thermal-averaged ghion-J/ip dissociation cross section as a 
function of the transverse momentum for Phard at (a) central and (b) forward rapidity 
regions. 



Without going into further details we simply quote the coupled differential equations that 
has to be solved to get the time dependence of ^ and phard [27] : 

^ :aA-'--^rno .S:^^-' (17) 



1 + ^ ^ r 27^(0 + 3(l + 07^'(0 

-^rPhard = ^7^ (4) , ^/i , ^\T,>^^\ i^^) 



1 + ePhard 27^(0 + 3(1 + OT^'iO 

where F = 2T(r)/(5f/) and f] = rj/s, rj is the shear viscosity co-efficient. In this model the 
time Tc has been calculated using the relation 7?.^/^(^) phard = Tc. We have used the same 
transverse profile for the hard momentum scale as in model I. A comparative study of the 
survival probability using the above described space-time models will be done. 



3 Results 

Let us first discuss the thermal averaged gluon dissociation cross section in anisotropic sys- 
tem. Eqs.(4)-(9) have been used for this purpose. The results are displayed in Fig. (1). 
Fig. (la) and Fig. (lb) correspond to Pt=0 and Pt=8 GeV respectively for a set of values of 
the anisotropy parameter. It is seen that the cross section decreases with ^ for Phard up to 
~ 500 MeV and then increases as compared to the isotropic case (^ = 0)(see in Fig.(la)). 
Similar feature has been observed in Fig. (lb) for higher P^ where the cross section starts to 
increase beyond Phard ~ 200 MeV. 

As mentioned before, the dissociation cross section depends on the direction of propaga- 
tion {6p) of the quarkonium with respect to the anisotropy axis. This dependence is shown in 
Fig. (2). We find marginal dependence in this case. These observations will have important 
consequences while calculating the survival probability (see later). 
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Figure 4: (Color online) Time evolutions of (a) the anisotropy parameter ^ and (b) the hard 
momentum scale Phard in the two space time models described in the text. 
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Figure 5: (Color online) The survival probability of J/V" in an anisotropic plasma at central 
rapidity, (a) corresponds to Tj=550 MeV, rj=0.7 fm/c and (b) is for Tj=440 MeV, rj=0.15 
fm/c. 
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In order to show the transverse momentum dependence for fixed hard momentum scale, 
we, in Fig. (3), present the dissociation cross section as a function of Pt of the J/ip for 
Phard = 300 McV and for two values of 6p. Again it is seen that the cross section first 
decreases with the anisotropy parameter upto Pt ~ 5(3) GeV for 9p = n/2{n/3) and we 
find larger increase away from the central rapidity region. This behavior might infiuence the 
survival probability which we consider next. 

Before calculating the survival probability let us examine the time evolutions of ^ and 
Phard in the space time models described earlier. This is needed to calculate the survival 
probability in an expanding plasma. The results are shown in Fig. (4). It is seen that the 
anisotropy parameter falls much rapidly compared to the case when model II is used (see 
Fig. (4a)). There is a narrow window in r where ^ dominates in case of model I. The coohng 
is slower in case of model II as can be seen from Fig. (4b). These observations have important 
consequence on the survival probability as we shall see. 

Eq.(lO) has been used to calculate the survival probability. For the space time model I we 
use Eq.(14) and for model II Eqs.(19) and (20) have been used for the time evolution of the 
anisotropy parameter ^ and the hard momentum scale Phard- Fig- (5) describes the survival 
probability for two different set (Set I and Set II) of initial conditions for a given direction of 
propagation of the J/ip with the anisotropy axis. For the Set 1 initial conditions, the results 
are same for the isotropic case and the two space time models used for the anisotropic media 
(see in Fig. (5a)). However, for the Set II initial conditions and for the same 9p, the results 
are different from each other as can be seen from Fig. (5b). More interestingly, we find an 
order of magnitude increase in the survival probability for the Set II initial conditions. This 
is because of the argument of the exponential in Eq.(lO). 

Next we consider the survival probability in forward rapidity region for the two sets of 
initial conditions and two space time models and compare it with that in the central rapidity 
region. Fig. (6a) ((6b)) shows the survival probability calculated using Set I (Set II) initial 
conditions. It is seen that for the Set I initial conditions, S{pt), in the forward rapidity 
region is marginally larger than the case when space time model II is used. It is seen from 
Fig. (6b) that for Set 11 initial conditions, S{pt) is marginally higher in the case of space time 
model I in the low px region at central rapidity. However, at forward rapidity, the survival 
probability is always larger in case of space time model I. 

In order to show the dependences of S{pt) on Tiso in the space time model I and on ri/s 
in space time model II we plot the survival probability for set I initial conditions in Fig. (7). 
It is observed that increasing Tiso lowers S{pt) while increasing r]/s enhances the survival 
probability. 

The results for the LHC energies for the two sets of space time models are shown in 
Fig. (8). In case of central rapidity region we do not find any difference in the results for the 
two space time models. But the results are marginally different in case of forward rapidity 
region. However, in the forward rapidity region we observe that the survival probability 
increases by a factor of 2 compared to the case of central rapidity for a given set of initial 
conditions. 
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Figure 6: (Color online) The survival probability of J/?/^ in an anisotropic plasma at central 
and forward rapidity regions corresponding to (a) Tj=550MeV, rj=0.7 fm/c and (b) Tj=440 
MeV, ri=0.15fm/c. 




Figure 7: (Color online) The survival probability of J j^l) in an anisotropic plasma for two 
different model at same initial temperature Tj = 500 MeV and initial time = 0.7 fm/c for 
different Tiso and 77/5. 
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Figure 8: (Color online) The survival probability of J/ip in an anisotropic plasma for (a) 
central and (b) forward rapidity region. The initial conditions are taken Tj = 820 Mev and 
Ti = 0.5 fm/c. 



4 Summary 

We have calculated gluon-J/ip dissociation cross section assuming pre-equilibrium momen- 
tum space anisotropy in the deconfined phase expected to be produced in relativistic heavy 
ion collisions. It is observed that the thermally weighted cross section is modified substan- 
tially in anisotropic plasma. To calculate the survival probability of the J/ip two sets of 
initial conditions and two different space models have been used both for RHIC and LHC 
energies. For set I initial conditions, in the central rapidity region we do not find any differ- 
ence in the survival probabilities calculated in the isotropic and anisotropic QGP at RHIC 
energies. However, changing the initial conditions, it is seen that the survival probability is 
lower in AQGP. Moreover, there is noticeable difference in the results obtained using differ- 
ent space time models. We also show that the results for the survival probability depend on 
the isotropization time (in model I) and rj/s (in model II). For the case of LHC with a given 
initial conditions the results marginally differ from each other in the forward rapidity region 
when two different space time models are used. It is also found that the results are not much 
sensitive to the direction of propagation of the J/ip with respect to the anisotropy axis. It is 
also demonstrated that the results are extremely sensitive to the initial conditions, in par- 
ticular, to the choice of the initial time. It is important to note that uncertainty may result 
from our assumption of chemical equilibrium. However, one can naively expect that finite 
chemical potentials should affect isotropic and anisotropic plasmas equally. So one expects 
that although the total yields could change, one would still see a sensitivity to the assumed 
isotropization or thermalization time. At leading order in the quark fugacity, the ratio of the 
isotropic to anisotropic result should be independent of the fugacity [28, 29]. We would also 
like to add that the consideration of transverse expansion may alter the results during the 
late stages of the collisions as has been observed in case of photon and dilepton transverse 
momentum distribution. However, transverse expansion is pronounced in the later stage 
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and its effect in tfie very early stage is minimal. Since in our case momentum anisotropy 
is an early stage phenomena, the effect will be negligible. We also not that treating the 
quarkonium binding energy as function of ^ and Phard might alter the present findings and 
this is worth investigating. 

It is to be noted that apart from this mechanism of suppression, there are other mecha- 
nisms by which J/ t/j can be suppressed [3] . All these possible processes should be taken into 
account and then be compared with the experimental data of transverse momentum distri- 
bution of J/ip to extract the isotropization time Tiso as has been done in case of photons [12] 
and nuclear modification of light hadrons [16]. 
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